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As an example for the interplay of structure, dynamics, and phase behavior of macromolecu- 
lar systems, this article focuses on the problem of bottle-brush polymers with either rigid or 
flexible backbones. On a polymer with chain length TV;,, side-chains with chain length TV are 
endgrafted with grafting density a. Due to the multitude of characteristic length scales and 
the size of these polymers (typically these cylindrical macromolecules contain of the order of 
10000 effective monomeric units) understanding of the structure is a challenge for experiment. 
But due to excessively large relaxation times (particularly under poor solvent conditions) such 
macromolecules also are a challenge for simulation studies. Simulation strategies to deal with 
this challenge, both using Monte Carlo and Molecular Dynamics Methods, will be briefly dis- 
cussed, and typical results will be used to illustrate the insight that can be gained. 



1 Introduction 

The so-called "bottle-brush polymers" consist of a long macromolecule serving as a "back- 
bone" on which many flexible side chains are densely grafted 1,2 . Since the chemical syn- 
thesis of such complex polymeric structures has become possible, these molecular bottle- 
brushes have found much interest for various possible applications: the structure reacts 
very sensitively to changes of solvent quality (due to change of temperature of the solu- 
tion, pH value, etc.), enabling the use of these molecules as sensors or actuators on the 
nanoscale 3 ' 4 . For some conditions these bottle-brush molecules behave like stiff cylindri- 
cal rods, and hence they can serve as building blocks of supramolecular aggregates, or 
show orientational order in solution as nematic liquid crystals do. On the other hand, these 
molecular bottle-brushes are also very "soft", i.e, they show only very small resistance to 
shear deformation; bottle-brush molecules of biological origin such as aggrecan which oc- 
curs in the cartilage of mammalian (including human!) joints are indeed held responsible 
for the excellent lubrication properties (reducing frictional forces) in such joints 5 . 

However, for being able to control the function of these complex macromolecules one 
must be able to control their structure, i.e., one must understand how the structure depends 
on various parameters of the problem: chain length Nf, of the "backbone", chain length N 
of the side chains, grafting density a of the side chains along the backbone, solvent quality, 
just to name the most important ones of these parameters. This is the reason why computer 




Figure 1. Snapshot pictures of simulated bottle-brush polymers with a flexible backbone containing JV;, = 131 
effective monomers. At each backbone monomer one side chain is grafted (grafting density cr = 1) and side 
chain lengths are N = 6 (a), TV = 12 (b), N = 24 (c) and N = 48 (d). The chains are described by the bond 
fluctuation model on the simple cubic lattice (see Sec. 2) and the only interaction considered is the excluded 
volume interaction between effective monomers. 

simulations are needed in order to understand such bottle-brushes: even in the case of 
good solvent, the structure of these objects is very complicated as the snapshot pictures 6,7 
(Fig. [TJ of the simulated models show, and there occur a multitude of characteristic length 
scales (Fig.O 6 ' 7 that are needed to describe the structure; it is very difficult to extract this 
information from experiment, and hence simulations are valuable here as they can yield a 
far more detailed picture. 

Of course, the large size of these bottle-brush polymers (even the coarse-grained mod- 
els described in the following sections contain of the order of several thousand or even 
more than 10000 effective monomers per polymer) is a serious obstacle for simulation, 
too; in fact, developing efficient models and methods for the simulation of macromolecu- 
lar systems is a longstanding and important problem 8 ' 9 . We hence shall address this issue 
of proper choice of both suitable models and efficient algorithms in the next sections. 

2 Bottle-brush polymers with rigid backbones in good solvents 

If one assumes the backbone of the bottle-brush to be completely stiff, the problem is re- 
duced to grafting side chains of length N to a straight line. Although this limiting case 
seems somewhat artificial, from the point of view that one wishes to model real systems, it 
is a very useful test case: there is no reason to assume that approximations that already fail 
in this rather simple limit become accurate for the more complicated case of flexible back- 
bones; moreover, this case is rather simple to simulate, and the analysis of the simulation 
data is relatively straightforward. 

For the study of this case, polymers were simply represented as standard self-avoiding 
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Figure 2. Multiple length scales are needed for the structural characterization of a bottle-brush polymer: a coarse- 
grained view describes this object as a flexible spherocy Under of length L cc and cross-sectional radius R CB , 
which is locally straight over a length t v ("persistence length"). A further length scale of interest is the end- 
to-end distance R e ^ of the backbone chain. This coarse-grained view is experimentally obtained by atomic 
force microscopy imaging of bottle-brush polymer adsorbed on substrates, for instance. A "microscope" with 
resolution on the scale of atoms would see the monomers of the backbone chains (connected by backbone bond 
vectors if,) and monomers of the side chains (connected by bond vectors £ s ). A mesoscopic length of interest 
then is the end-to-end distance R e of the side chains. 



walks on the simple-cubic lattice (i.e., beads of the chains are occupied lattice sites, con- 
nected by nearest-neighbor links on the lattice, and multiple occupancy of sites is forbid- 
den 8 ' 9 ). 

As a simulation method, pruned-enriched Rosenbluth methods (PERM) 10, 11 were 
used 12 . In the Rosenbluth method, all side chains are grown simultaneously step by step, 
choosing only from sites which are not yet occupied, and the statistical weight of the poly- 
mer configuration is computed recursively. In the PERM algorithm, one does not grow a 
single polymer at a time, but a large "population" of equivalent chains is grown simultane- 
ously but from time to time configurations with very low statistical weight are killed, and 
configurations with large statistical weight are "cloned" ("go with the winners"strategy 10 ). 
The advantage is that (unlike dynamic Monte Carlo algorithms 8 ' 9 ) this method does not 
suffer from "critical slowing down" when N gets large, and data for all N (up to the max- 
imum value studied) are obtained in a single simulation (see 11 for more details). 

Now one hypothesis popular in most experimental studies (e.g. 13 ' 14 ) is the factorization 
approximation 15 for the structure factor S(q) (that describes the small angle scattering 
intensity of neutrons, X-rays or light from dilute solutions of bottle-brushes polymers) into 
a contribution due to the backbone (S(q)) and due to the side chains (S s (q)) 



S(q) « S b (q)S s (q) » S b (q)S cs (q) (D 
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Figure 3. a) Log-log plot of the normalized scattering function of the backbone St(q)/Li,, being the back- 
bone length, if, = Ni,a where a(= 1) is the lattice spacing. Three choices of L;, are shown, b) Log-log plot of 
the scattering from all the side chain monomers, for N = 50, and two choices of the grafting density c) Log-log 
plot of the scattering due to interference contributions from monomers in side chains and the backbone (for the 
same parameters as in b). d) Log-log plot of qS cs (q) vs. q in the range 0.01 < q < 1, comparing the actual 
scattering due to the side chains with the prediction resulting from Eq. (2j> using p(r) as actually recorded in the 
simulation 12 . 



where q is the absolute value of the scattering wave vector. In the last step the side chain 
scattering was approximated by S cs (q). This cross sectional scattering is in turn approx- 
imated related to the radial density p cs (r) perpendicular to the cylinder axis at which the 
side chains are grafted as (the constant c ensures proper normalization 12 ) 



S cs (q) = 




Obviously, Eq. fl} neglects interference effects due to correlations between the monomer 
positions in the side chains and in the backbone. Writing the scattering due to the 
monomers in the side chains as the average of a square {Eq. (O} ignores correlations 
in the occupation probability in the z-direction along the bottle-brush axis. All such corre- 
lations do contribute to the actual structure factor, of course, when it is computed from its 
definition, 

TVtot Mot 

S ^ = aT~ E X>ft) c fo)> Sin(q\n - fjD/iqlfi - fj\) . (3) 

^ l= l 3 =1 

Here c(ri) is an occupation variable, c(fi) = 1 if the site i is occupied by a bead, and 
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qR g 

Figure 4. Structure factor S(q) of bottle-brush polymer obtained 13 from scattering experiment mapped to the 
simulated model (cf. text) by requiring that the total gyration radius R g is matched, to fix the translation factor 
for the length scale 7 . 
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Figure 5. Log-log plot of the density profile p C s(r) vs. r, for the systems shown in Fig. [4] For fitting the 
experimental data, a Gaussian form for p C s(r) was assumed 13 . 



zero otherwise, and the sums extend over all A/tot monomers (in the side chains and in the 
backbone). 

In the simulation, one can go beyond experiment by not only recording the total scat- 
tering S(q), but also the contributions Sb{q), S s (q) and S cs (q) individually (Fig.0. One 
can see that the approximation S(q) ~ Sb(q)S s (q) leads to a relative error of the order 
of a few % only (cf. the difference in ordinate scales between parts b) and c), while the 
assumption S s (q) w S cs (q) [with Eq. ©J leads to appreciable errors at large q (see part 
d)). Therefore the use of Eqs. ((T), © to analyze experiments can lead to appreciable errors, 
when one wishes to predict 12 the radial density distribution p(r). 
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Figure 6. Schematic diagram of states of a bottle-brush polymer with a rigid backbone under poor solvent condi- 
tions in the plane of variables x = jV 1 / 2 ^ and y = TV 1 / 2 ^ — T/8), as proposed by Scheiko et al. 19 , The lines 
indicate (smooth) crossovers between the different states 19 , which we have characterized by snapshots from our 
simulations 20 . For further explanations see text. 



3 Bottle-brushes with flexible backbones in good solvents 



In the case of flexible backbones, the PERM algorithm no longer provides an efficient 
sampling of configuration space; there "dynamic" Monte Carlo algorithms are used, but 
with "unphysical moves" allowing chain intersections and large configurational changes in 
a single step, as achieved by the Pivot algorithm 8 ' 9 . Using the bond fluctuation model 8 9 , 
the "L26 algorithm" 16 was used 17 for local moves allowing bond intersection and never- 
theless respecting excluded volume. In this way well-equilibrated data for systems with up 
to Nf, — 259 monomers at the backbone and up to N = 48 monomers in the side chains 
could be studied (for a = 1). Examples of snapshot pictures have already been given in 
Fig. ID 

Adjusting the physical meaning of the lattice spacing to correspond to 0.263 nm, the 
structure factor of the simulated model matches almost perfectly an experimental result 13 
(for slightly different numbers of chemical monomers, N™ p — 400, iV cxp = 62; this dif- 
ference is irrelevant, since there is no one-to-one correspondence between covalent bonds 
and the "effective bonds" of the model) see 7 Fig. |4] From the simulation one can directly 
extract the cross-sectional density profile p cs (r) and compare it 7 to the approximate exper- 
imental result (Fig. [5]l, which was obtained from fitting S(q) to Eqs. ([T), 13 - One sees 
that the analysis of the experiment could predict roughly correctly the distance on which 
the profile p cs (r) decays to zero, but does not account for its precise functional form. A 
further interesting finding 7, 17 is the result that the persistence length l v (cf. Fig5]i depends 
strongly 7 on Nf, (at least if one uses the textbook definitions 18 ), and hence is not a useful 
measure of local chain stiffness. 
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Figure 7. Log-log plot of wavelength A (dimensionless by normalization with cr) versus grafting density a (di- 
mensionless by normalization with the radius of gyration of the side chains, Rg). The wavelength was extracted 
either from the z-dependence of a normalized correlation function at radial distance r ss 3 Lennard- Jones units 
a (stars) or extracted from the peak position of its Fourier transform S(q) (squares) [inset]. All data refer to 
temperature T = 1.5 (in Lennard Jones units). Blue and red symbols refer to N = 35, orange and green ones to 
N = 50. Inset shows S(q) vs. q for N = 35 and grafting densities a > 0.38. 



4 Bottle-brushes in poor solvents 



While isolated single chains in dilute solution collapse to dense globules 18 when the tem- 
perature T is lowered below the Theta temperature T — 9, for a bottle-brush the constraint 
of the chains being grafted to a backbone leads to an interesting diagram of states, consid- 
ering 1 — T/9 and a as variables (Pig. [6j> 19 . Near T = 9 and oN 1 ! 2 — > the side chains 
collapse to dense globules when N 1/2 (l - T/9) > 1. For T » 9 and N 1/2 a > 1 one 
has a cylindrical structure (as in the previous sections) which collapses to a dense cylinder 
when iV 1 / 2 (l — T/9) 1. However, interestingly at intermediate grafting densities and 
N 1 / 2 (l-T/9) > 1 a laterally inhomo geneous "pearl necklace-structure" was predicted . 

This problem can neither be studied efficiently by the PERM method nor by the bond 
fluctuation algorithm of Sec. 3 - both methods get very inefficient for dense polymer con- 
figurations. Thus, instead Molecular Dynamics simulations of the standard Grest-Kremer- 
type bead-spring model 8,9 were carried out 20 . Fig. [7] shows 20 characteristic wavelength A 
plotted vs. grafting density a. One finds that for small a a trivial periodicity given by I /a 
occurs, while for aR g > 0.2 (R g = (R 2 ) 1 / 2 ) the periodicity is essentially independent of 
cr, given by about that A w 8R g . Of course, mapping out the phase behavior precisely for 
a wide range of N, a and T is a formidable problem, due to very long relaxation times of 
the collapsed states in Fig. [6] 
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5 Concluding Remarks 

The above examples have illustrated that polymers with complex architecture pose chal- 
lenging problems with respect to their structure and thermodynamics. Investigation of their 
relaxation behavior is another interesting problem, to be tackled in the future. The above 
examples have also illustrated the need to carefully adjust both the model and the algorithm 
to the problems one wants to deal with. 
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